Evolution of T cells in the cancer-resistant naked mole-rat.

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αßT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

Urocortin-3 neurons in the mouse perifornical area promote infant-directed neglect and aggression.

While recent studies have uncovered dedicated neural pathways mediating the positive control of parenting, the regulation of infant-directed aggression and how it relates to adult-adult aggression is poorly understood. Here we show that urocortin-3 (Ucn3)-expressing neurons in the hypothalamic perifornical area (PeFAUcn3) are activated during infant-directed attacks in males and females, but not other behaviors. Functional manipulations of PeFAUcn3 neurons demonstrate the role of this population in the negative control of parenting in both sexes. PeFAUcn3 neurons receive input from areas associated with vomeronasal sensing, stress, and parenting, and send projections to hypothalamic and limbic areas. Optogenetic activation of PeFAUcn3 axon terminals in these regions triggers various aspects of infant-directed agonistic responses, such as neglect, repulsion, and aggression. Thus, PeFAUcn3 neurons emerge as a dedicated circuit component controlling infant-directed neglect and aggression, providing a new framework to understand the positive and negative regulation of parenting in health and disease.

The Tug1 lncRNA locus is essential for male fertility.

BACKGROUND: Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function. RESULTS: We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function. CONCLUSIONS: Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Transitions in paternal social status predict patterns of offspring growth and metabolic transcription.

One type of parental effect occurs when changes in parental phenotype or environment trigger changes to offspring phenotype. Such nongenetic parental effects can be precisely triggered in response to an environmental cue in time-locked fashion, or in other cases, persist for multiple generations after the cue has been removed, suggesting multiple timescales of action. For parental effects to serve as reliable signals of current environmental conditions, they should be reversible, such that when cues change, offspring phenotypes change in accordance. Social hierarchy is a prevalent feature of the environment, and current parental social status could signal the environment in which offspring will be born. Here, we sought to address parental effects of social status and their timescale of action in mice. We show that territorial competition in seminatural environments affects offspring growth. Although dominant males are not heavier than nondominant or control males, they produce faster growing offspring, particularly sons. The timing, effect-size, and sex-specificity of this association are modulated by maternal social experience. We show that a change in paternal social status is sufficient to modulate offspring weight: from one breeding cycle to the next, status-ascending males produce heavier sons than before, and status-descending males produce lighter sons than before. Current paternal status is also highly predictive of liver transcription in sons, including molecular pathways controlling oxidative phosphorylation and iron metabolism. These results are consistent with a parental effect of social experience, although alternative explanations are considered. In summary, changes in paternal social status are associated with changes in offspring growth and metabolism.

Single-cell transcriptomics of the naked mole-rat reveals unexpected features of mammalian immunity.

The immune system comprises a complex network of specialized cells that protects against infection, eliminates cancerous cells, and regulates tissue repair, thus serving a critical role in homeostasis, health span, and life span. The subterranean-dwelling naked mole-rat (NM-R; Heterocephalus glaber) exhibits prolonged life span relative to its body size, is unusually cancer resistant, and manifests few physiological or molecular changes with advancing age. We therefore hypothesized that the immune system of NM-Rs evolved unique features that confer enhanced cancer immunosurveillance and prevent the age-associated decline in homeostasis. Using single-cell RNA-sequencing (scRNA-seq) we mapped the immune system of the NM-R and compared it to that of the short-lived, cancer-prone mouse. In contrast to the mouse, we find that the NM-R immune system is characterized by a high myeloid-to-lymphoid cell ratio that includes a novel, lipopolysaccharide (LPS)-responsive, granulocyte cell subset. Surprisingly, we also find that NM-Rs lack canonical natural killer (NK) cells. Our comparative genomics analyses support this finding, showing that the NM-R genome lacks an expanded gene family that controls NK cell function in several other species. Furthermore, we reconstructed the evolutionary history that likely led to this genomic state. The NM-R thus challenges our current understanding of mammalian immunity, favoring an atypical, myeloid-biased mode of innate immunosurveillance, which may contribute to its remarkable health span.

Murine single-cell RNA-seq reveals cell-identity- and tissue-specific trajectories of aging.

Aging is a pleiotropic process affecting many aspects of mammalian physiology. Mammals are composed of distinct cell type identities and tissue environments, but the influence of these cell identities and environments on the trajectory of aging in individual cells remains unclear. Here, we performed single-cell RNA-seq on >50,000 individual cells across three tissues in young and old mice to allow for direct comparison of aging phenotypes across cell types. We found transcriptional features of aging common across many cell types, as well as features of aging unique to each type. Leveraging matrix factorization and optimal transport methods, we found that both cell identities and tissue environments exert influence on the trajectory and magnitude of aging, with cell identity influence predominating. These results suggest that aging manifests with unique directionality and magnitude across the diverse cell identities in mammals.

eIF2B activator prevents neurological defects caused by a chronic integrated stress response.

The integrated stress response (ISR) attenuates the rate of protein synthesis while inducing expression of stress proteins in cells. Various insults activate kinases that phosphorylate the GTPase eIF2 leading to inhibition of its exchange factor eIF2B. Vanishing White Matter (VWM) is a neurological disease caused by eIF2B mutations that, like phosphorylated eIF2, reduce its activity. We show that introduction of a human VWM mutation into mice leads to persistent ISR induction in the central nervous system. ISR activation precedes myelin loss and development of motor deficits. Remarkably, long-term treatment with a small molecule eIF2B activator, 2BAct, prevents all measures of pathology and normalizes the transcriptome and proteome of VWM mice. 2BAct stimulates the remaining activity of mutant eIF2B complex in vivo, abrogating the maladaptive stress response. Thus, 2BAct-like molecules may provide a promising therapeutic approach for VWM and provide relief from chronic ISR induction in a variety of disease contexts.

Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region.

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

Correction to: A window into extreme longevity; the circulating metabolomic signature of the naked mole-rat, a mammal that shows negligible senescence.

The original version of this article unfortunately contained an error.

A window into extreme longevity; the circulating metabolomic signature of the naked mole-rat, a mammal that shows negligible senescence.

Mouse-sized naked mole-rats (Heterocephalus glaber), unlike other mammals, do not conform to Gompertzian laws of age-related mortality; adults show no age-related change in mortality risk. Moreover, we observe negligible hallmarks of aging with well-maintained physiological and molecular functions, commonly altered with age in other species. We questioned whether naked mole-rats, living an order of magnitude longer than laboratory mice, exhibit different plasma metabolite profiles, which could then highlight novel mechanisms or targets involved in disease and longevity. Using a comprehensive, unbiased metabolomics screen, we observe striking inter-species differences in amino acid, peptide, and lipid metabolites. Low circulating levels of specific amino acids, particularly those linked to the methionine pathway, resemble those observed during the fasting period at late torpor in hibernating ground squirrels and those seen in longer-lived methionine-restricted rats. These data also concur with metabolome reports on long-lived mutant mice, including the Ames dwarf mice and calorically restricted mice, as well as fruit flies, and even show similarities to circulating metabolite differences observed in young human adults when compared to older humans. During evolution, some of these beneficial nutrient/stress response pathways may have been positively selected in the naked mole-rat. These observations suggest that interventions that modify the aging metabolomic profile to a more youthful one may enable people to lead healthier and longer lives.

Functional circuit architecture underlying parental behaviour.

Parenting is essential for the survival and wellbeing of mammalian offspring. However, we lack a circuit-level understanding of how distinct components of this behaviour are coordinated. Here we investigate how galanin-expressing neurons in the medial preoptic area (MPOAGal) of the hypothalamus coordinate motor, motivational, hormonal and social aspects of parenting in mice. These neurons integrate inputs from a large number of brain areas and the activation of these inputs depends on the animal's sex and reproductive state. Subsets of MPOAGal neurons form discrete pools that are defined by their projection sites. While the MPOAGal population is active during all episodes of parental behaviour, individual pools are tuned to characteristic aspects of parenting. Optogenetic manipulation of MPOAGal projections mirrors this specificity, affecting discrete parenting components. This functional organization, reminiscent of the control of motor sequences by pools of spinal cord neurons, provides a new model for how discrete elements of a social behaviour are generated at the circuit level.

Benign and tumor parenchyma metabolomic profiles affect compensatory renal growth in renal cell carcinoma surgical patients.

BACKGROUND AND OBJECTIVES: Pre-operative kidney volume is an independent predictor of glomerular filtration rate in renal cell carcinoma patients. Compensatory renal growth (CRG) can ensue prior to nephrectomy in parallel to tumor growth and benign parenchyma loss. We aimed to test whether renal metabolite abundances significantly associate with CRG, suggesting a causative relationship. DESIGN, SETTING, PARTICIPANTS, AND MEASUREMENTS: Tissue metabolomics data from 49 patients, with a median age of 60 years, were previously collected and the pre-operative fold-change of their contra to ipsi-lateral benign kidney volume served as a surrogate for their CRG. Contra-lateral kidney volume fold-change within a 3.3 +/- 2.1 years follow-up interval was used as a surrogate for long-term CRG. Using a multivariable statistical model, we identified metabolites whose abundances significantly associate with CRG. RESULTS: Our analysis found 13 metabolites in the benign (e.g. L-urobilin, Variable Influence in Projection, VIP, score = 3.02, adjusted p = 0.017) and 163 metabolites in the malignant (e.g. 3-indoxyl-sulfate, VIP score = 1.3, adjusted p = 0.044) tissues that significantly associate with CRG. Benign/tumor fold change in metabolite abundances revealed three additional metabolites with that significantly positively associate with CRG (e.g. p-cresol sulfate, VIP score = 2.945, adjusted p = 0.033). At the pathway level, we show that fatty-acid oxidation is highly enriched with metabolites whose benign tissue abundances strongly positively associate with CRG, both pre-operatively and long term, whereas in the tumor tissue significant enrichment of dipeptides and benzoate (positive association), glycolysis/gluconeogenesis, lysolipid and nucleotide sugar pentose (negative associations) sub-pathways, were observed. CONCLUSION: These data suggest that specific biological processes in the benign as well as in the tumor parenchyma strongly influence compensatory renal growth.

The Multipartite Mitochondrial Genome of Enteromyxum leei (Myxozoa): Eight Fast-Evolving Megacircles.

Myxozoans are a large group of poorly characterized cnidarian parasites. To gain further insight into their evolution, we sequenced the mitochondrial (mt) genome of Enteromyxum leei and reevaluate the mt genome structure of Kudoa iwatai. Although the typical animal mt genome is a compact, 13-25 kb, circular chromosome, the mt genome of E. leei was found to be fragmented into eight circular chromosomes of ∼23 kb, making it the largest described animal mt genome. Each chromosome was found to harbor a large noncoding region (∼15 kb), nearly identical between chromosomes. The protein coding genes show an unusually high rate of sequence evolution and possess little similarity to their cnidarian homologs. Only five protein coding genes could be identified and no tRNA genes. Surprisingly, the mt genome of K. iwatai was also found to be composed of two chromosomes. These observations confirm the remarkable plasticity of myxozoan mt genomes.

New Perspectives on Genomic Imprinting, an Essential and Multifaceted Mode of Epigenetic Control in the Developing and Adult Brain.

Mammalian evolution entailed multiple innovations in gene regulation, including the emergence of genomic imprinting, an epigenetic regulation leading to the preferential expression of a gene from its maternal or paternal allele. Genomic imprinting is highly prevalent in the brain, yet, until recently, its central roles in neural processes have not been fully appreciated. Here, we provide a comprehensive survey of adult and developmental brain functions influenced by imprinted genes, from neural development and wiring to synaptic function and plasticity, energy balance, social behaviors, emotions, and cognition. We further review the widespread identification of parental biases alongside monoallelic expression in brain tissues, discuss their potential roles in dosage regulation of key neural pathways, and suggest possible mechanisms underlying the dynamic regulation of imprinting in the brain. This review should help provide a better understanding of the significance of genomic imprinting in the normal and pathological brain of mammals including humans.

Regulation of alternative splicing at the single-cell level.

Alternative splicing is a key cellular mechanism for generating distinct isoforms, whose relative abundances regulate critical cellular processes. It is therefore essential that inclusion levels of alternative exons be tightly regulated. However, how the precision of inclusion levels among individual cells is governed is poorly understood. Using single-cell gene expression, we show that the precision of inclusion levels of alternative exons is determined by the degree of evolutionary conservation at their flanking intronic regions. Moreover, the inclusion levels of alternative exons, as well as the expression levels of the transcripts harboring them, also contribute to this precision. We further show that alternative exons whose inclusion levels are considerably changed during stem cell differentiation are also subject to this regulation. Our results imply that alternative splicing is coordinately regulated to achieve accuracy in relative isoform abundances and that such accuracy may be important in determining cell fate.

Genomic insights into the evolutionary origin of Myxozoa within Cnidaria.

The Myxozoa comprise over 2,000 species of microscopic obligate parasites that use both invertebrate and vertebrate hosts as part of their life cycle. Although the evolutionary origin of myxozoans has been elusive, a close relationship with cnidarians, a group that includes corals, sea anemones, jellyfish, and hydroids, is supported by some phylogenetic studies and the observation that the distinctive myxozoan structure, the polar capsule, is remarkably similar to the stinging structures (nematocysts) in cnidarians. To gain insight into the extreme evolutionary transition from a free-living cnidarian to a microscopic endoparasite, we analyzed genomic and transcriptomic assemblies from two distantly related myxozoan species, Kudoa iwatai and Myxobolus cerebralis, and compared these to the transcriptome and genome of the less reduced cnidarian parasite, Polypodium hydriforme. A phylogenomic analysis, using for the first time to our knowledge, a taxonomic sampling that represents the breadth of myxozoan diversity, including four newly generated myxozoan assemblies, confirms that myxozoans are cnidarians and are a sister taxon to P. hydriforme. Estimations of genome size reveal that myxozoans have one of the smallest reported animal genomes. Gene enrichment analyses show depletion of expressed genes in categories related to development, cell differentiation, and cell-cell communication. In addition, a search for candidate genes indicates that myxozoans lack key elements of signaling pathways and transcriptional factors important for multicellular development. Our results suggest that the degeneration of the myxozoan body plan from a free-living cnidarian to a microscopic parasitic cnidarian was accompanied by extreme reduction in genome size and gene content.

Quantitative and functional interrogation of parent-of-origin allelic expression biases in the brain.

The maternal and paternal genomes play different roles in mammalian brains as a result of genomic imprinting, an epigenetic regulation leading to differential expression of the parental alleles of some genes. Here we investigate genomic imprinting in the cerebellum using a newly developed Bayesian statistical model that provides unprecedented transcript-level resolution. We uncover 160 imprinted transcripts, including 41 novel and independently validated imprinted genes. Strikingly, many genes exhibit parentally biased--rather than monoallelic--expression, with different magnitudes according to age, organ, and brain region. Developmental changes in parental bias and overall gene expression are strongly correlated, suggesting combined roles in regulating gene dosage. Finally, brain-specific deletion of the paternal, but not maternal, allele of the paternally-biased Bcl-x, (Bcl2l1) results in loss of specific neuron types, supporting the functional significance of parental biases. These findings reveal the remarkable complexity of genomic imprinting, with important implications for understanding the normal and diseased brain.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes.

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.

The operonic location of auto-transcriptional repressors is highly conserved in bacteria.

Bacterial genes are commonly encoded in clusters, known as operons, which share transcriptional regulatory control and often encode functionally related proteins that take part in certain biological pathways. Operons that are coregulated are known to colocalize in the genome, suggesting that their spatial organization is under selection for efficient expression regulation. However, the internal order of genes within operons is believed to be poorly conserved, and hence expression requirements are claimed to be too weak to oppose gene rearrangements. In light of these opposing views, we set out to investigate whether the internal location of the regulatory genes within operons is under selection. Our analysis shows that transcription factors (TFs) are preferentially encoded as either first or last in their operons, in the two diverged model bacteria Escherichia coli and Bacillus subtilis. In a higher resolution, we find that TFs that repress transcription of the operon in which they are encoded (autorepressors), contribute most of this signal by specific preference of the first operon position. We show that this trend is strikingly conserved throughout highly diverged bacterial phyla. Moreover, these autorepressors regulate operons that carry out highly diverse biological functions. We propose a model according to which autorepressors are selected to be located first in their operons in order to optimize transcription regulation. Specifically, the first operon position helps autorepressors to minimize leaky transcription of the operon structural genes, thus minimizing energy waste. Our analysis provides statistically robust evidence for a paradigm of bacterial autorepressor preferential operonic location. Corroborated with our suggested model, an additional layer of operon expression control that is common throughout the bacterial domain is revealed.

Evolutionary models accounting for layers of selection in protein-coding genes and their impact on the inference of positive selection.

The selective forces acting on a protein-coding gene are commonly inferred using evolutionary codon models by contrasting the rate of nonsynonymous substitutions to the rate of synonymous substitutions. These models usually assume that the synonymous substitution rate, Ks, is homogenous across all sites, which is justified if synonymous sites are free from selection. However, a growing body of evidence indicates that the DNA and RNA levels of protein-coding genes are subject to varying degrees of selective constraints due to various biological functions encoded at these levels. In this paper, we develop evolutionary models that account for these layers of selection by allowing for both among-site variability of substitution rates at the DNA/RNA level (which leads to Ks variability among protein-coding sites) and among-site variability of substitution rates at the protein level (Ka variability). These models are constructed so that positive selection is either allowed or not. This enables statistical testing of positive selection when variability at the DNA/RNA substitution rate is accounted for. Using this methodology, we show that variability of the baseline DNA/RNA substitution rate is a widespread phenomenon in coding sequence data of mammalian genomes, most likely reflecting varying degrees of selection at the DNA and RNA levels. Additionally, we use simulations to examine the impact that accounting for the variability of the baseline DNA/RNA substitution rate has on the inference of positive selection. Our results show that ignoring this variability results in a high rate of erroneous positive-selection inference. Our newly developed model, which accounts for this variability, does not suffer from this problem and hence provides a likelihood framework for the inference of positive selection on a background of variability in the baseline DNA/RNA substitution rate.